Low-energy relativistic effects and nonlocality in time-dependent tunneling 
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We consider exact time-dependent analytic solutions to the Schrodinger equation for tunneling 
in one dimension with cut off wave initial conditions at t = 0. We obtain that as soon as t 7^ 
the transmitted probability density at any arbitrary distance rises instantaneously with time in a 
linear manner. Using a simple model we find that the above nonlocal effect of the time-dependent 
solution is suppressed by consideration of low-energy relativistic effects. Hence at a distance xo 
from the potential the probability density rises after a time to = xo/c restoring Einstein causality. 
This implies that the tunneling time of a particle can never be zero. 

PACS numbers: 03.65,73.40. Gk 
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Recent technological achievements as the possibility of 
constructing artificial quantum structures at nanomet- 
ric scales [EJ or manipulating individual atoms [0 have 
stimulated a great deal of work at both an applied and 
fundamental level. In particular, studies on tunneling 
have addressed, among others, the controversial question 
of the traversal time of a particle through a classically 
forbidden region Q . The above considerations have mo- 
tivated a renewed attention to the time-dependent treat- 
ments of quantum tunneling. From the theoretical side, 
most of these works are based on the numerical analy- 
sis of the time-dependent Schrodinger equation with the 
initial condition of a Gaussian wave packet Q] . A com- 
mon feature in most of these approaches is that the initial 
wave packet extends through all space. As a consequence 
the initial state, although it is manipulated to reduce as 
much as possible its value along the tunneling and trans- 
mitted regions, contaminates from the beginning the tun- 
neling process. In the literature, however, one also finds 
a number of approaches to time-dependent tunneling, pi- 
onnered by Stevens [^) , that in fact circumvent the above 
situation using a cut off wave as initial state |^|,|[0J§1 ■ 

Our approach is a generalization to an arbitrary poten- 
tial H of the Moshinsky shutter ||. Moshinsky consid- 
ered the solution of the time-dependent free Schrodinger 
equation with the initial condition, at t = 0, of a plane 
wave of momentum k confined in the half-space region 
x < to the left of a perfectly absorbing shutter lo- 
cated at x — 0. The sudden opening of the shutter at 



time t = 0, allows the plane wave solution to propagate 
freely along the region x > |l(J. Moshinsky showed 
that as the time t goes to infinity, the solution to the 
problem tends to the stationary solution. He also found 
that both the current and the probability density for a 
fixed value of the distance Xq as a function of t, present 
oscillations near the wavefront, situated at to = xq/v. 
He named this phenomenon diffraction in time, in anal- 
ogy to the well known phenomenon of optical diffraction. 
Recently an observation of diffraction in time has been 
reported ]0]. If we put a potential barrier in the re- 
gion < x < L with the same initial condition as above, 
then we may have a convenient model to analyze tun- 
neling times by measuring at what time the probability 
density rises from zero. However, as pointed out by Hol- 
land [12 for the free case, and by Garcfa-Calderon and 
Rubio ||], for the case of a potential, the solution of the 
time-dependent Schrodinger equation for a cut off initial 
plane wave has a nonlocal character. This means that 
if initially there is a zero probability for the particle to 
be at x > 0, as soon as t ^ 0, there is instantaneously 
a finite, though very small, probability to find it at any 
point x > 0. This implies a zero tunneling time for some 
particles. 

In this work we address the issue of the behaviour of 
the time-dependent solution to the Schrodinger equation 
for tunneling through a potential barrier using a cut off 
wave as initial condition. Our aim is to analyze the non- 
local behaviour of the time-dependent transmitted solu- 
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tion at early times. We also study low-energy relativistic 
effects by solving the Klein-Gordon equation for a model 
potential. The implication of our findings for the tunnel- 
ing time problem is briefly discussed. 

For the sake of simplicity in our approach we con- 
sider the instantaneous removal of the shutter. This may 
be seen as a kind of 'sudden approximation' to a shut- 
ter opening with finite velocity, where the treatment is 
more involved even for the free case E^]. In a recent 
paper we have shown that the transmitted solution for 
the Schrodinger case for tunneling through an arbitrary 
potential barrier may be written as a free term solution 
plus a infinite sum resonance transient terms associated 
with the S-matrix poles of the problem [||. This cor- 
responds to solve the time-dependent Schrodinger equa- 
tion for a potential V(x) that vanishes outside the region 
< x < L, with the initial condition, 



ip s (x, k,t = 0) = 



e ikx , x < 
0, x > 0. 



(1) 



The transmitted solution for the region x > L reads 

oo 

^ s (x,k,t)=T(k)M(x,k,t)-iJ2 T n M (x,k n ,t), (2) 

n 

where T{k) stands for the transmission amplitude of the 
problem, T n = u n (0)u n (L)exp(—ik n L)/(k — k n ), is given 
in terms of the resonant eigenfunctions u n (x) and com- 
plex S-matrix poles k n ; and the Moshinsky functions 
M(x, k, t) and M(x, k n ,t) are denned as, 



M(x,q,t) = -e^ imx2 / 2ht h v2 evtc(y), 
where the argument y is given by 
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In the above two equations q stands either for k or k n . 
In the absence of a potential the solution given by Eq. 
(^J) becomes the solution for the free case obtained by 
Moshinsky §J, 



ip°(x,k,t) = M(x,k,t). 



(5) 



As discussed by Moshinsky, the initial condition given by 
Eq. (Q) refers to a shutter that acts as a perfect absorber 
(no reflected wave) . One can also envisage a shutter that 
acts as a perfect reflector. In such a case the initial wave 
may be written as, 



i/} s (x, k,t = 0) = 



-ikx 



, x < 



(6) 



0, x > 0. 

The transmitted solution for the region x > L now reads, 



tp s (x,k,t) =T(k)M(x,k,t) -T(-k)M(x,-k,t) 

OO 

-2ikJ2T n M{x,k„,t), (7) 

n 

where T n = u n (0)u n (L)exp(~ik n L) / (k 2 — k 2 ,) . The solu- 
tion for the free case with the reflecting initial condition 
is 



ip°(x,k,t) = M(x,k,t) - M(x,-k,t). 



(8) 



The exact solutions given by Eqs. (g) and (Q), corre- 
sponding, respectively, to absorbing and reflecting initial 
cut off waves, involve each a contribution proportional 
to the free case solution and then an infinite sum in- 
volving the S-matrix poles, {fc„}, and resonant states, 
{u n (x)}, of the system. As shown in ref. ||, at very 
long times, the terms M{x,k ni t) that appear in the 
above equations vanish. The same occurs for M(x, —k, t) 
while, as shown firstly in ref. Q, M(x,k,t) tends to 
the stationary solution. Hence, at long times, each of 
the above exact solutions go into the stationary solution 
T(k)exp[i(kx ~ Et/h)}. 

At very short times, for a given x > L, the argument 
of M(x, k,t), given by Eq. (Q) with q — k, becomes very 
large and in fact becomes independent of the value of 
k, y w exp(— iTr/4:)[m/2ht] 1 ^ 2 x. Since for very large y, 
M(y) ~l/y it follows that M(x, k,t) goes like t 1 / 2 . 
As discussed also in ref. Bj, the functions dependent on 
the poles, M(x, k n ,t), behave in a similar fashion pro- 
vided the value oft — to is sufficiently small to guarantee, 
for a fixed x = L, that L 3> h\k n \t/m. Since the distribu- 
tion of the complex S-matrix poles on the k-planc fulfills 
[ | fci| < | A;2 1 - - - < \k n \..., one sees that as t becomes smaller 
and smaller there will be more and more values of n for 
which the corresponding M functions goe like t 1 / 2 as do 
all the rest of M functions associated with smaller values 
of n. In the appropriate limit as t — > and n — > oo, the 
corresponding M function then vanishes as t 1 / 2 . Conse- 
quently for x > L, the solutions given by Eqs. (0) and 
(0) are proportional to t 1 / 2 , namely, 



A 



ip s (x,k,t) i 1/2 , (x>L) 



(9) 



where A a constant. Note that at t = the solution 
vanishes in accordance with the initial condition. It is 
not difficult to see that Eq. (||) will hold also for a cut 
off initial condition that is something between the ini- 
tial conditions considered above, and more generally, for 
a wave packet formed by a linear combination of cut 
off waves. Eq. (||) tell us that the probability density 
at any distance x from the potential will rise instanta- 
neously with time. This intriguing nonlocal behaviour 
implies that an ideal detector will measure a zero tunnel- 
ing time. The existence of action-at-a-distance effects in 
the time-dependent Schrodinger equation should not in 
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principle pose any conceptual difficulties since the treat- 
ment is non-relativistic. However one could ask whether 
the above nonlocal behaviour arises because the initial 
condition is a cut off wave. In order to answer the above 
question we consider low-energy relativistic effects by 
solving the Klein-Gordon equation with a cut off wave as 
initial condition for a simple potential model. Moshinsky 
H solved the Klein- Gordon equation for the free shut- 
ter problem with the initial condition of a cut off plane 
wave in the region x < and showed that the probabil- 
ity density at a point x > is nonzero only after a time 
to > xq/c, with c the velocity of light. To our knowledge 
a numerical analysis of this solution has not yet been per- 
formed. We would like to learn also how the relativistic 
solution is affected at early times by tunneling through 
a potential. 

A potential that has been widely used in studies on 
time-dependent tunneling is the square barrier, char- 
acterized by a height Vq and a width L. This poten- 
tial has an infinite set of S-matrix poles situated at in- 
creasing energies on top of the barrier. There is, how- 
ever, a simpler potential model that is more amenable 
for a relativistic treatment. This is the delta potential 
V(x) = b s 8(x). The solution corresponding to the time- 
dependent Schrodinger equation has been obtained by El- 
berfeld and Kleber using a delta- function propagator jl4| . 
One can also follow a derivation by Laplace transforming 
directly the time-dependent Schrodinger equation of the 
problem using the initial condition given by Eq. (|l|) Jl5[ . 
Defining p 2 = 2ims/h the Laplace transformed solution 
ip s (x, s) for the region x > reads, 



ip s (x,k,s) 



irn 



h (p + ib) (p — k) ' 



(10) 



where b = mb s /h 2 . After a simple partial fractions 
decomposition the inverse Laplace transform yields for 
x > 0, 

ip s s (x,k,t) = T(k)M(x,k,t) + R(k)M(x,-ib;t). (11) 

where T(k) and R(k) stand for the transmission and re- 
flection amplitudes for the stationary situation, T(k) = 
kj(k + ib) and R(k) = ib/(k + ib). Note that here in- 
stead of an infinite number of S-matrix poles the only 
S-matrix pole corresponds to an antibound state located 
at k a = —ib. At a very short times one can easily see 
that ipg(x, k,t) goes like t 1 / 2 fulfilling also, as the square 
barrier, Eq. (JsJ) . 

The shutter problem for the Klein-Gordon equation 
with the delta potential V(x) = b r 5(x) requires the solu- 
tion of 



dx 2 



~^m 2 ' 

-[b r 8(x) + fjL 2 ]ip s r (x, k r , t) 



■^- 1 ripf.(x,k r ,t) = ——ipf.(x,k r ,t) 



(12) 



iftr (%, k r ,t = 0) = ' 



ik r x 



x < 
0, x > 



(13) 



where we define E 2 = k 2 +[i 2 and k r = fc(l — (fc/^) 2 )~ 1//2 . 
Note that E r is given in units of the reciprocal length, 
i.e., E r = E/hc. The condition given by Eq. jl^ ) 
follows from the fact that for t < 0, when the shut- 
ter was closed, we had on the left side of the shutter, 
i//(x,k r ,t) = exp[i(fc r a; — E r ct)], for x < and a vanish- 
ing value for x > 0. By direct application of the Laplace 
transform method one gets a set of differential equations 
corresponding to the regions x > and x < 0. In order to 
derive an expression for the transmitted wave function, 
we have to consider the matching conditions to take into 
account the discontinuity of the wave function derivatives 
at x — 0, obtaining the Laplace-transformed solution, 



(s — icE r ) 



-qx/c 



2 (q + b c) (q + ik r c) 



(14) 



where q — [s 2 + ^c 2 ] 1 / 2 and bo = b r /2. Using the 
Bromwich contour to evaluate the inverse Laplace trans- 
form of Eq. (H) it is convenient to make the change 
of variable || — iu = (q + s)/(/ic). In this form q = 
ificiu -1 — u)/2 and as a consequence the branch points 
at s = ±iy,c go into an essential singularity at u = and 
two simple poles located on the lower half of the complex 
it-plane. After separating into partial fractions one then 
may evaluate the resulting integrals by standard complex 
variable techniques to obtain the wave function, 



i> s r (x,k r ,t) 



BD[ip°(x,- 
0, 



■BCi>° r {x,-ib ,t)- 
ibo,t)]*, t > x/c 
t < x/c 



(15) 



where A — k r /(k r + ibo), B = ibo/(k r + ibo), C — 
(e + E r )/(2e), and D = (e - E r )/(2e), and also e = 
(/i 2 - bl) 1 ' 2 . In Eq. @, the function i()°(x,k r ,t) is 
the solution of the free Klein-Gordon case |)jl6|, namely, 



e i(k rX -E r ct) + l j (^)_ 

j OO 

^ (X, k r ,t) = { E E/«f MV), t > X/C (16) 
n=0 

0, t < x/c 

where J„ (rj) stands for the Bessel function of order n and, 



ct 



ct 



1/2 



, 77 = n(c 2 t 2 



E r ). 

(17) 



where /1 = mc/h. with the initial condition given by, 



The expressions, ip^(x,—ibo,t) and (ip® (x, — ibo, t))* in 
Eq. dla) have the same form as the free solution in 
Eq. ( |l6| ) with k r replaced by —ibo- Asymptotically for 
very long times in the solution ^(x, k r ,t), given by Eq. 
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(|l5|), the terms ip^{x, —ib ,t) and (ip®(x,—ibo,t))* van- 
ish, while the term i/}®(x,k r ,t) goes into the stationary 
solution exp[i(fc r a; — E r ct)]. 

To exemplify the above results Fig. [j] exhibits the very 
short time behaviour of the probability density for the 
delta potential at a fixed distance x = L = 0.3A . One 
sees that the Schrodinger description (broken line), ob- 
tained from Eq. ( pd| ) with parameters b s = 2.0 eV — A 
and E = 0.01 eV, yields an instantaneous response with 
time while the relativistic solution, calculated using Eq. 
(|i"5|), starts after to == L/c. This tell us something rel- 
evant: The nonlocal behaviour of the Schrodinger de- 
scription is due to its non-relativistic nature. The non- 
local behavior of the Schrodinger solution would result 
from the fact that in a non-relativistic description there 
is no restriction on the velocity of some components of 
the initially confined wave function. The sharp rela- 
tivistic wavefront of height 0.25 in Fig. [j] follows as a 
consequence of the initial condition given by Eq. (13). 
This jump occurs also in the free case and may be ob- 
tained analytically J9| . For an initial function of the type 
exjp(ik r x) + exp(za)exp(— ik r x), (x < 0), with a an ar- 
bitrary phase, the peak height will be a function of a. 
In particular for a reflecting initial condition, (a = tt), 
the solution starts smoothly from zero at to — L/c. It 
might be of interest to mention that in fully relativistic 
quantum field theories Hegerfeldt |l7j has pointed out 
that the sudden opening of a shutter may lead to vio- 
lation of Einstein causality, i.e., no propagation faster 
than light. This author has argued that the difficulty is 
of a theoretical nature and has discussed some ways to 
solve it. Our relativistic model satisfies Einstein causal- 
ity. The inset to Fig. |l] shows that at longer times the 
above two solutions approach each other, both presenting 
the characteristic transient behaviour near the 'classical' 
wavefront at x = vt, which in our example occurs at a 
very short time. Our analysis has a consequence of inter- 
est for the tunneling time problem. Since the probability 
density rises with time after a time to — xo/c, it implies 
that the tunneling time of a particle can never be zero, 
contrary to some claims in the literature || . 

Thus we can see that a proper description of the quan- 
tum mechanical propagation for the transmitted solu- 
tion, even at low energies, strictly requires a relativis- 
tic treatment. However, since the corresponding solu- 
tions are practically identical up to the relativistic cut 
off, at t — L/c, suggests that the Schrodinger description 
is quite accurate provided the velocity components larger 
than c are omitted. 
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FIG. 1. Plot of \ipi{x,t)\ 2 (dashed line) and 
\il>r{x,t)\ 2 (solid line), respectively, for the Schrodinger and 
Klein-Gordon solutions for a delta potential, as a function 
of time at a fixed distance at early times and at long times 
(inset). See text. 
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